Magnetism in graphene induced by single-atom defects 
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We study from first principles the magnetism in graphene induced by single carbon atom defects. 
For two types of defects considered in our study, the hydrogen chemisorption defect and the vacancy 
defect, the itinerant magnetism due to the defect-induced extended states has been observed. Cal- 
culated magnetic moments are equal to 1 \xb per hydrogen chemisorption defect and 1.12—1.53 \ib 
per vacancy defect depending on the defect concentration. The coupling between the magnetic mo- 
ments is either ferromagnetic or antiferromagnetic, depending on whether the defects correspond to 
the same or to different hexagonal sublattices of the graphene lattice, respectively. The relevance 
of itinerant magnetism in graphene to the high- Tc magnetic ordering is discussed. 

PACS numbers: 61.72.Ji, 75.75.+a, 81.05.Uw 



I. INTRODUCTION 

The last two decades were marked with the dis- 
coveries of new allotropic modifications of carbon and 
related nanostructures. Graphene, the single two- 
dimensional sheet of graphite, is the starting point for 
many carbon nanomaterials, which are commonly called 
nanographites. These materials, being diverse in atomic 
structure, display a wide range of electronic properties. 
Magnetism of carbon materials 1 is of particular inter- 
est since in current technological applications magnetic 
materials are based on d and / elements. New carbon- 
based magnetic materials would greatly extend the limits 
of technologies relying on magnetism. Even more promis- 
ing is the application of such materials in the design of 
nanoscale magnetic and spin electronics devices. 

While ideal graphite and carbon nanotubes are in itself 
nonmagnetic, experimental observations of magnetic or- 
dering are often explained by the presence of impurities 2 , 
boundaries 3 -^ or defects 5 -^^. Defects in nanographites 9 
can be created intentionally by irradiating material with 
electrons or ion o 10 ' 11 ' 12 ' 13 . By manipulating the condi- 
tions of irradiation it is possible to tune, in a flexible way, 
the properties of the carbon-based material s 14 ' 15 : 16 ' 17 . 
Examples of simple defects in nanographites are single 
atom vacancies and hydrogen chemisorption defects. The 
former defect type is produced upon the irradiation with 
high energy particles 6 * while the latter is the major out- 
come of the hydrogen plasma treatment^. The com- 
mon feature of both types of defects is that only one 
carbon atom is removed from the ir conjugation network 
of the graphene sheet. The single-atom defects on the 
graphene lattice give rise to quasilocalized states at the 
Fermi leve l 10 ' 18 ' 19 ' 20 . The graphene lattice is the bipar- 
tite lattice. It can be viewed as two interpenetrating 
hexagonal sublattices of carbon atoms (labeled a and 
(3). When a defect is created in the a lattice, only the p z 
orbitals of carbon atoms in the (3 sublattice contribute 
to the quasilocalized state, and vice versa. These states 
extend over several nanometers around the defects form- 



ing characteristic (v^x ^3)^30° superstructures recog- 
nized in STM images. Analyzing the position and the 
orientation of the superstructures one can precisely lo- 
cate the defect and determine the sublattice to which it 
belong s 18 ' 21 . The fact that quasilocalized states lie at the 
Fermi level suggests that itinerant (Stoner) magnetism 
can be induced by the electron exchange instability. It 
has been argued recently, that the Stoner ferromagnetism 
with high Curie temperatures, Tc, can be expected for sp 
electron systems^ 2 . On the other hand, the Ruderman- 
Kittel-Kasuya-Yoshida (RKKY) coupling^ 3 - of localized 
magnetic moments in graphene is too weak to result 
in high Tq*^ Another recent work has shown that the 
RKKY-type interactions in graphene are suppressed. 25 

In this work we investigate using first principles ap- 
proaches the itinerant magnetism originating from the 
quasilocalized states induced by single-atom point de- 
fects in graphene. The results obtained can eventually 
be extended, with some precautions, to defects in other 
nanographites. 



II. COMPUTATIONAL METHODS AND 
MODELS 

The model system consists of a periodic two- 
dimensional superlattice of defects in graphene (Fig. QJ,). 
The supercell size can be varied resulting in different dis- 
tances d between the neighbor defects on the superlattice 
and, thus, in different defect concentrations. The results 
can be further extrapolated to the cases of low defect 
concentrations. For the chosen supercell, the resulting 
distance between neighbor defects is about 3na cc where 
a cc =1.42 A is the carbon-carbon distance in graphene. 
The corresponding number of carbon atoms per unit cell 
is 6n 2 . Our investigation is restricted to the cases with 
ri=2— 6. The largest system considered (n=6) is char- 
acterized by about 25 A separation between neighbor 
defects, which corresponds to a defect concentration of 
0.5%. 



2 




FIG. 1: (a) Definition of the extendable two-dimensional 
hexagonal lattice of defects in the graphene sheet. The unit 
cell and the Voronoi cell are shown as full and dashed lines, 
respectively. The defective atom is labeled by the triangle ac- 
cording to the orientation of the defect state y/3xy/3 super- 
structure. The size of the supercell shown here corresponds to 
9a cc separation between neighbor defects (n=3). (b) Struc- 
ture of the hydrogen chemisorption defect, (c) Structure of 
the vacancy defect, (d) Hexagonal closest packing (n=l) of 
vacancy defects with the corresponding unit cell. 
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FIG. 2: (a) Calculated magnetic moment per defect vs. the 
separation d between the neighbor hydrogen chemisorption 
defects (•) and the neighbor vacancy defects (o). (b) Fermi 
energies («,o), majority spin (A, A) and minority spin (T,V) 
bands maxima versus the defect separation d for the hydrogen 
chemisorption defects (filled symbols) and the vacancy defects 
(open symbols). The Fermi level of ideal graphene is taken as 



Density functional theory calculations were performed 
using the SIESTA code i 26 i 27 The generalized gradient ap- 
proximation exchange-correlation density functional of 
Perdew, Burke and Ernzerhof (PBE)2£ was employed. 
All calculations were performed in the spin-unrestricted 
manner using the diagonalization-based method for solv- 
ing Kohn-Sham equations. The shifted Monkhorst-Pack 
grids 2 ^ corresponding to a cutoff of 100 Bohr were used 
to sample the Brillouin zone in two dimensions. Atomic 
positions and cell dimensions were relaxed. The numer- 
ical atomic orbital basis set of single-^ plus one polar- 
ization function (SZP) quality was used for the whole 
range of models studied. All calculations for the mod- 
els with n=2— 4 were reproduced using the basis set of 
double-^ plus one polarization function (DZP) quality. 
For all electronic structure quantities discussed in this 
study (magnetic moment, Fermi levels and band max- 
ima), there is a good agreement between the results of 
the two basis sets, despites the slight overestinration of 
the C— C bond length found in the SZP calculations. 
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FIG. 3: Density of states plots for the systems with the hy- 
drogen chemisorption defects (a) and with the vacancy defects 
(b) (n = 4). The dashed line shows the density of states of 
the ideal graphene. Labels indicate the character of the defect 
states. 



III. DISCUSSION OF RESULTS 

In the following we present our results for the two types 
of defects mentioned above. The structure of the hydro- 
gen chemisorption defect is shown in Fig.[]J>. This defect 
is characterized by the slight protrusion of the hydro- 
genated carbon atom and the very small displacement 
of all other neighbor carbon atoms^^ The single atom 
vacancy defect in graphene is nearly planar (Fig. [TJ;). 
The local three-fold symmetry breaks down due to the 
Jahn- Teller distortion induced by the reconstruction of 
two dangling bonds left after removing the carbon atom. 



This gives rise to the in-plane displacement of other car- 
bon atoms in the graphene lattice^^i The third dan- 
gling bond is left unsaturated providing a contribution 
of magnitude 1 /is to the intrinsic magnetic moment 
of the defect. For the case of the vacancy type defect 
(Fig.[TJi) in the closest packing geometry (n=l) no single 
six-membered ring remains. This interesting structure 
can be considered as yet another hypothetical allotropic 
modification of carbon for which one may expect a high 
specific magnetic moment. 

Magnetism induced by the presence of quasilocalized 



defect states ipd(r) has been observed in the case of 
both defect types. The hydrogen chemisorption defect 
gives rise to the strong Stoner ferromagnetism^ with 
a magnetic moment of 1 fiB per defect at all studied 
concentrations (Fig. [2k). The flat defect bands give 
rise to the very narrow peaks (W<0.2 eV) which are 
necessary for the stability of magnetic ordering at high 
temperatures^ The defect band maxima for the major- 
ity spin and the minority spin components lie, respec- 
tively, lower and higher than the the Fermi levels for both 
defective and ideal graphene (Fig. [5Jd). The hydrogen 
chemisorption motif is charge neutral and spin-polarized 
in the wide range of defect concentrations. On the 
contrary, fractional magnetic moments and weak Stoner 
ferromagnetism^ have been observed for the vacancy- 
type defect models. A magnetic moment of 1.15 \xb has 
been predicted for the closest packing of vacancy type 
defects (n—1) (Fig. [TJi), while for smaller defect con- 
centrations the magnetic moment was found to vary in 
the range of 1.45—1.53 \xb per defect (Fig. [5k). For the 
vacancy- type defect, the total magnetic moment is de- 
termined by the contribution (1 hb) of the localized sp 2 
dangling bond state (atom 1 in Fig. [3} and the contribu- 
tion (<1 (j,b) of the extended defect state, ipd(,r) (labeled 
Pz in Fig. GO). 

The width of the defect state bands and the overall 
modification of the band structure are larger in the case 
of the vacancy type defects (Fig. [5b). The partial spin 
polarization of ipd(r) (filled majority spin band and half- 
filled minority spin band) is explained by the self-doping 
(charge transfer from the bulk to ipd{r)), which arise from 
the stabilization of the defect state. The stabilization of 
vacancy defect extended states is possible in the case of a 
significant coupling between the second nearest neighbor 
atoms belonging to the same sublattice^ In the case of 
the vacancy defect, the indirect coupling is justified by 
the formation of the covalent bond between the two car- 
bon atoms V (Fig. [4)3) that follows the defect reconstruc- 
tion. No such bond is possible in the case of hydrogen 
chemisorption. Thus, the character of the defect-induced 
magnetism depends on the possibility of covalent bond- 
ing between the second nearest neighbor atoms due to the 
reconstruction. This provides an interesting opportunity 
for tailoring magnetic properties of materials. 

The defect state exchange splitting, defined as the dif- 
ference between the corresponding majority spin and mi- 
nority spin band maxima, decreases as the defect con- 
centration decreases. This is not surprising since the 
degree of the localization of the defect states depends 
on the defect concentration^ At the lowest studied de- 
fect concentration of 0.5%, the exchange splitting de x 
were found to be 0.23 eV and 0.14 eV for the hydro- 
gen chemisorption and vacancy defects, respectively. In 
the latter case, the splitting is smaller due to the partial 
spin polarization of the defect band. Since in both cases 
de x >k B T for T «300 K, the Stoner theory predicts T c 
above room temperature for defect concentrations of the 
order of 1%. The decrease of the Stoner theory Tc due to 
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FIG. 4: (Color online) Spin density projection (in /^s/a.u. 2 ) 
on the graphene plane around (a) the hydrogen chemisorption 
defect (A) and (b) the vacancy defect in the a sublattice. 
Carbon atoms corresponding to the a sublattice (o) and to 
the /3 sublattice (•) are distinguished. Simulated STM images 
of the defects are shown in (c) and (d), respectively. 



spin wave excitations is expected to be ineffective for the 
case of carbon based materials^ At low concentrations 
the magnetism in defective nanographitcs is expected to 
be sensitive to the variations of the Fermi energy result- 
ing from self-doping, presence of other defects or applied 
bias, and to the disorder-induced broadening. 

The distributions of the electron spin magnetization 
density in the vicinity of both types of defects clearly 
show the characteristic \/3x\/3 patterns also observed 
for the charge density in the STM experiments. For the 
hydrogen chemisorption defect the projection of the spin 
density (Fig. 0k) on the graphene plane clearly shows 
three-fold symmetry. For the vacancy type defect the 
symmetry is broken due to the Jahn- Teller distortion 
(Fig. |4p) . The localized magnetic moment associated to 
the dangling bond of atom 1 can also be observed. The 
simulated STM images (Fig. [4J; and (Hi) based on our 
calculations agree with experimental observations i 10 ' 18 i 19 
The distribution of the electron spin density is repre- 
sented in Fig. [5] (n—6 model) by means of the Mulliken 
spin populations averaged over ith nearest neighbors to 
the defect atom. The spin populations show a damped 
oscillation behavior as a function of the nearest neighbor 
index and, therefore, of the distance to the defect. The 
magnetization pattern is explained by the fact that the 
defect state is distributed over the sites of the sublattice 
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FIG. 5: Dependence of the spin populations averaged over ith 
nearest neighbors of the hydrogen chemisorption defects (•) 
and the vacancy defects (o). The spin population for the first 
nearest neighbor atoms of the vacancy defect (0.39) are out 
of scale due to the contribution of the localized sp 2 dangling 
bond state. 



complementary to the one in which the defect was cre- 
ated (i. e. over the odd nearest neighbors), and shows 
a power law decay^ The major positive contribution to 
the electron spin density is defined by the exchange split- 
ting of the defect states. In addition, the exchange spin- 
polarization effect (i.e. the response of the fully popu- 
lated valence bands to the magnetization of the defect 
states) results in a negative spin density on the even 
nearest neighbor sites and in the enhancement of a pos- 
itive spin density on the odd nearest neighbor sites (see 
Fig. QJ. A similar phenomenon takes place in the case 
of the neutral bond length alternation defect states (neu- 
tral solitons) in one-dimensional polyene chains i 33 i The 
calculated magnitude of the negative spin-polarization is 
«l/3 of the positive spin populations on the neighbor 
sites in the vicinity of the defect site. This is close to the 
ratio observed for the trans-polyacetylene^ 3 - The magni- 
tudes of the spin populations are lower in the case of the 
vacancy defect because of the fractional spin-polarization 
of the defect band. 

According to the Stoner picture, the magnetic order- 
ing is driven by the exchange energy E x ~ — ^\ M? 
with Mj being the magnetization of the p z orbital of the 
jth carbon atom^ Ferromagnetic ordering is the only 
possibility for the magnetism originating from quasilo- 
calized states induced by defects in the same sublat- 
tice because of the non-oscillating behavior of both (i) 
Mj within the same sublattice and (ii) the indirect 
(RKKY) coupling due to the semimetallic properties 
of graphene^ On the contrary, for the case of defect 
states in different sublattices, E x is minimized when the 
coupling is antiferromagnetic. In this case, the mech- 
anism of the exchange coupling is defined by the indi- 
rect spin-polarization effect. The strength of the cou- 
pling between the defect-induced magnetic moments lo- 
cated in different sublattices depends on the defect con- 




FIG. 6: (Color online) Spin density distribution in the sys- 
tem with three hydrogen chemisorption defects (two defects 
in sublattice a (>) and one in sublattice /3 (<)). The same 
scale as in Fig. [4] applies. 

centration since E x ~ — J2 i M? . The contribution of 
the magnetic moment associated with a single defect is 
~ E 3 \Mrj)\ 4 ~ log- 2 (N), where 1/N is the 
defect concentration.- To further illustrate this point, 
we calculated the ground state magnetic configuration 
of the system with three close hydrogen chemisorption 
defects using the DFT approach. We found that in the 
ground state configuration two defects in the a sublat- 
tice are coupled ferromagnetically with each other and 
antiferromagnetically with the third defect in the j3 sub- 
lattice (Fig. [6|). The resulting magnetic moment of this 
system amounts to 2—1=1 fiB in agreement with the sec- 
ond Lieb theorem^ and the characteristic V3x\/3 pat- 
terns of the magnetization density associated with indi- 
vidual defects can be recognized. In nanographite ma- 
terials with defects present with an equal probability in 
both sublattices, the overall correlation of the magnetic 
moments is expected to be antiferromagnetic. The an- 
tiferromagnetic ordering was experimentally observed in 
carbon nanohornsi 36 ' 37 

In addition, we found that the proton spin of 
chemisorption hydrogen atom shows a strong magnetic 
coupling to the electron in the quasilocalized defect state. 
Even for the lowest defect concentration studied we find 
the magnitude of Fermi contact hyperfine interaction of 
«35 Gauss. One can consider the nuclear spin quantum 
computer implementation originally proposed by Kane 3 - 
based on carbon nanostructures. 



IV. CONCLUSIONS 

In conclusion, our calculations reveal the itinerant 
magnetism triggered by simple defects in graphene and 
stable over the wide range of concentrations. It is no- 
table, that the itinerant magnetism does not require the 
presence of highly reactive unsaturated dangling bonds. 
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Both ferromagnetic and antiferromagnetic scenarios of 
the magnetic correlation are possible with the second be- 
ing more probable for truly disordered systems. The re- 
construction of vacancy defects was found to be responsi- 
ble for the partial suppression of magnetic moments and 
for the broadening of defect bands. We argue that the 
defect-induced itinerant magnetism is responsible for the 
experimentally observed high-Tc ferromagnetism of irra- 
diated graphite samples. 
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